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ABSTRACT 

Using the new evolutionary code Pegase.3, we made an evolutionary spectral synthesis of 
the optical-IR-submm spectral energy distribution of two distant (z = 3.8) radio galaxies, 
4C 41.17 and TN J2007-1316. These two radio galaxies were selected from the HeRGE 
(Herschel Radio Galaxies Evolution) Project in particular for their faint AGN contribution 
and because they show evidence of a large stellar contribution to their bolometric luminosity. 
Pegase.3 coherently models the reprocessing of the stellar luminosity to dust emission allow- 
ing us to build UV-to-IR-submm spectral energy distribution libraries which then can be used 
to fit spectral energy distributions in the observer's frame. Our principal conclusion is that a 
single stellar population is insufficient to fit the spectral energy distributions of either radio 
galaxy. Our best fits are a sum of two evolving stellar populations - a recent starburst plus an 
old population - plus the thermal emission from an active galactic nucleus (which provides a 
good fit to the mid-IR emission). The two stellar components are: i) a massive 10 11 M0) 
starburst ^30 Myrs after formation, which is required simultaneously to fit the far-IR Her- 
schel to submmdata and the optical dataii) an older massive 10 n ~ 12 M0) early-type galaxy 
population, 1.0 Gyr old, which is required principally to fit the mid-IR Spitzer/IRAC data. 
A young population alone is insufficient because an evolved giant-star population produces a 
1 /im rest-frame peak which is observed in the IRAC photometry. This discovery confirms that 
much of the stellar populations in high redshift radio galaxies are formed by massive starbursts 
in the early universe. Gas-rich mergers and/or jet-cloud interactions are favored for triggering 
the intense star formation necessary to explain the properties of their spectral energy distri- 
butions. The discovery of similar characteristics in two distant radio galaxies suggests that 
multiple stellar populations, one old and one young, may be a generic feature of the luminous 
infrared radio galaxy population as a whole. 

Key words: stars:formation-galaxies: evolution,active,stellar content -cosmology: 
miscellaneous-infrared: galaxies -radio continuum: galaxies 



1 INTRODUCTION 

Large scale cosmological simulations offer the prospect of 
constraining theories of galaxy formation and probing the appro- 
priateness of cosmological parameters. However, hydrodynamical 
and N-body simulations meet various difficulties in fitting obser- 
vations. The reasons for these difficulties is that there are many 
physical processes and timescales which are not well understood 
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such as: the appropriate initial conditions for the models, the 
detailed structure of dark matter halos, the timescale over which 
mass is accumulated into the halo and onto the galaxy, how stars 
form in detail, the relative role of major and minor mergers and 
their importance over cosmic time, the effect of the density of 
the galaxian environment and the impact of stellar and active 
galactic nuclei (AGN) feedback. Radio galaxies are known to 
be massive elliptical ty p es which formed in t h e early Universe 
jvan Breugel et al. 19981 : IPentericci et al. 200 it Izirm et al. 20031 : 
lLacv et al. 201 lh . The puzzling Hubble K - z diagram, which 
is bounded on the bright side principally by powerful radio 
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galaxies jLillv and Longair 1984 |Pe Breuck et al. 2002ah . 
allowed us to identify in the comoving rest-frame the 
population of elliptical galaxies. This population is young 
but already massive (10 12 M B ) out to high z =4 redshifts 
dRocca-Volmerange et al. 2 004). On the other hand the mid-IR, 
far-IR and sub-millimeter emission from distant radio galaxies 
were interpreted as sta rb ursts mainly initiat ed by major mergers 



dPe Breuck et al 20ld: llvison et al. 20081; llvison et al. 20121 ; 



lEngel et al. 201ol ; ISevmour et al. 201(1 ISevmour et al. 2012b . 
These dichotomous results: the old population from the Hubble 
K - z diagram and the young population from the need for high 
rates of star formation in radio galaxies imply that we are required 
to simultaneously follow the passive evolution of the galaxy as 
well as that of the on-going starburst to explain the overall spectral 
energy distribution (SED). To do this simultaneously requires a re- 
fined spectrophotometric model such as our new version of Pegase. 
The code Pegase. 3 offers the possibility to interpret observations 
in the observer's frame by using UV-to-submm synthetic libraries 
built for a wide range of ages. The template SEDs take into account 
the local z=0 templates corrected for evolution of the stellar pop- 
ulation (e-correction) and redshift (k-correction). In the present 
paper we tackle a precise, heretofore still debated key issue-namely 
whether the most distant galaxies, including radio galaxies, are 
pure starbursts initiated by on-going mergers or are forming stars 
continuously on timescales longer than that generally appropriate 
for mergers. In radio galaxies, there is also the additional difficulty 
of the thermal emission from the AGN (perhaps a torus) and of the 
scattered light from the AGN which both contribute to the overall 
spectral ener gy distribution, esp ecially in the mid-IR and UV 
respectively dDrouart et al. 2012b . We chose two radio galaxies in 
particular, 4C 41.17 and TN J2007-1316, both at z=3.8, because 
of their small AGN contribution to the overall SED. We apply the 
following two selection criteria: (i) relatively faint scattered light 
contribution in the rest-frame UV as determined from polarimetric 
measurements. Polarimetric measurements are the only way to 
ensure that the continuum is mostly of stellar origin. Indeed 
those two galaxies are the only cases where a c lear signature of 
photospheric absorption lines has been detected dPev et al. 19971) 
(ii) faint contribution of the AGN in the rest-frame 8-12/um range 
of the SED. To fit the SED, one and two stellar components are 
successively tried from an automatic procedure which minimizes 
the x 2 - We focus on the main properties (type, mass, age and 
star-formation history) of the two stellar populations. In § 2, we 
present the observations of the two radio galaxies for which we 
have built the continuous optical-IR-submm SEDs after correcting 
for instrumental effects (differing apertures, flux calibrations, etc.). 
Pegase. 3 and the fitting method are described in § 3. § 4 presents 
the best fit results of SEDs in the observer's frame obtained using 
two evolved stellar populations plus a simple AGN model. Finally, 
we provide some thoughts on the implications of these results in § 5 
and our conclusions in the last section. The adopted cosmological 
parameters are Ho = lOkm.s^' Mpc^ 1 ,Q. M = 0.3, Qa = 0.7. 



2 OBSERVATIONS 

The two radio galaxies 4C 41.17 (z = 3.80) and TN J2007-1316 
(z — 3.84) are part of the sample of ultra-steep spectrum radio 
sources (USS; a < -1.3; S v <x v") designed specifically to in- 
crea se the probability of d iscovering distant, z > 3, radio galax- 
ies dPe Breuck et al. 2000h . The two targets were specifically se- 
lected because of their relatively low level of AGN activity in the 



rest-frame UV/optical and their strong photospheric and dust sig- 
natures respectively seen in their optical continuum emission and 
in the cold dust-grain dominated far-IR emission. The Spitzer High 
Redshift Radio Galaxies (SHzRG) sample, which these two galax- 
ies are also part of, were augmented wit h A"-band photometry 
dSevmour et al. 20071 ; ISevmour et al. 20 08). Most of the comple- 
mentary optical and submm photometr y and data reduction p ro- 
cedures have been described previously toe Breuck et al. 20101) . In 
addition, both galaxies have recently been observed by Herschel 
with the PACS (100 and 160 pm) an d SPIRE (250, 350 and 500 pm) 
instruments, see also for 4C 41.17 dWvlezalek et al. 20121) . and are 
part of the HeRGE (Herschel Radio Galaxy Evolution) Project 
dSevmour et al. 20T2I) . Optical to submm spectral energy distribu- 
tions, sampling the rest-frame photospheric (UV/optical/near-IR) 
stellar emission to emission from cold grains (far-IR) are carefully 
built paying specially attention to aperture effects, flux calibration, 
and flux transmission of the filters used for the observations and 
removing the contribution to the fluxes from strong emission lines. 
For the two galaxies the synchrotron emission is considered as neg- 
ligible in the wavelength domain hereafter analyzed due to the steep 
radio spectra of the cores. 



2.1 The radio galaxy 4C41.17 

The radio galaxy 4C 41. 17 (z = 3.800 ± 0.00 3) is the archetype 
of distant radio galaxies (Cham bers et al. 19901) . A detailed multi- 
frequency radio analysis with the VLA and MERLIN shows that 
the radio structur e is associated with high resolutio n optical imag- 
ing HST data dCarilli et al. 19941 ; iMilev et al. 19921) . The galaxy is 
detected with a good signal-to-noise ratio and is spatially resolved 
at the 0.1" HST resolution. Based on both the Cj//i to Cr y line 
ratio and the strength of the C/y line, dBicknell et al. 200d) sug- 
gested that there is an interaction of the high-powered jet with a 
dense cloud in the halo of 4C 41.17. Such an interaction through 
high speed shocks (~ 1000 km/s), leads to strong compression 
of some of the gas, inducing star formation. The bright, spatially 
extended rest-frame UV continuum emission from this galaxy, 
aligned with the radio axis, is unpolarised (f 20- < 2.4%) and shows 
tell-tale stellar absorpti on features indicative of on-going star for- 
mation dDev et al. 199% . Two massive (M dy „ = 6 x 10 10 M Q ) com- 
ponents have been identified from CO(J=4-3) observations with 
the IRAM interferometer a nd interpreted as an evidence of merg- 
ers dPe Breuck et al. 2 005). However, the CO (J=l-0) emission 
line is not detected and the lower limits obtained are typical of 
molecular gas in starbursts JPapadopoulo s et al. 2005b . To com- 
plete the broad wavelength coverage of the SED of 4C 41.17, we 
use the total of the flux from 4C 41.17 within 6 arcsec field in the 
/f j-band, which is basically free o f emission from strong optical 
emission lines 1 Grah am et al. 19941) . For the submillimeter emis- 
sion w e use the results f rom the JCMT dPunlop et al. 1994)b and 
IRAM dChini et al. 19941) . 

The velocity and dispersion fields w ere mapped with t he in- 
tegral field spectrograph TIGER/CFHT dAdam et al. 1997h with 
a 0.61 arcsec spatial sampling. Radial velocities are essentially 
negative as the velocity field shows bow shocks with high ve- 
locity dispersions. Perpendicular to the main radio axis, a South- 
west extension of velocity = -600 kms~' and velocity dispersion 
= 2000 kms~' is visible. Table Q] contains the photometry data we 
used to construct the SED of 4C 41.17. The upper limits and data 
not used by the fitting procedure are also provided for completeness 
and are also show in subsequent figures. 
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Table 1. Photometric data, in red on Fig. [5] and Fig. [4] of the radio 
galaxy 4C41.17 including new Herschel/SPIKE and PACS observations 
(o-(J2000) = 06 h 50 m 52.098 s , <5(J2000) = +41°30'30.53"). -tot- in column 
4 means that galaxy photometry is totally measured within the instrument 
aperture. Complementary data, lower rows, in green on Fig. [f]and Fig.[4] are 
not used for fits, but are provided for completeness. 
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2.2 The radio galaxy TN J2007-1316, z=3.84 

TN J2007-1316 (or WN J2007-1316)is also an ultra steep spec- 
trum radio source s elected from the 365 MH z TEXAS surve y 
JPouglas et al. 199fjh and the NVSS survey l lCondon et al. 1998h . 
It was also detect ed in the 352 MHz WIS H survey and is part of 
the USS sample jDe Breuck et al. 2002bl) . Spectropolarimetry of 
TN J2007-1316 obtained with FORS1 on the VLT on 3 nights 
from 2002 May 6 to 8(Fig. Q}. For these observations, the 300V 
grism was used with a 1" wide slit oriented North-South and had 
a resolution of ~ 10 A (FWHM). Conditions were photometric 
with seeing between 0.5-0.9 arcsecs. During each night, four equal- 
length exposures were taken using four different orientations of the 
half-wave plate (0°, 45°, 22°5, 67°5). The total integration time was 
27600s. The initial data reduction followed standard procedures in 
the NOAO IRAF package. Because the object is relatively faint in 
the single 35 or 40 minute exposures (each split into o- and e-rays), 
we needed to pay particular attention to extracting the same phys- 
ical apertures in each of the individual spectra. We therefore first 
contructed the total 2-dimensional o- and e-ray spectra, and defined 
1x2.4 arcsec 2 wide apertures from these spectra. We then used the 
aperture and trace of these spectra to extract the 12 o-ray and 12 e- 
ray spectra, using the same linear dispersion (2.645 A /pix) in order 
to calculate the polarization in equally-sized spectral bins. Finally, 
we combined these spectra using the median to construct the eight 
spectra needed to calculate the pola rization vector. We followed 
the procedures of dVernet et al. 20011) to calculate the polarization 
percentage and position angle. We checked the polarisation angle 
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Figure 1. VLT spectropolarimetry of TN J2007-1316. the top panel is di- 
vided in three parts. The upper shows the total intensity spectrum, the mid- 
dle plot is the percentage polarization, and the bottom panel the polarization 
angle, as function of observed wavelengths. Continuum bins are denoted by 
a dot. Vertical error bars denote 1 <x uncertainties. The grey histograms show 
the data with a 150 A binning. The bottom panel is the 2D spectrogram, 
showing evidences of stellar signatures. 



offset between the half-wave plate coordinate and sky coordinates 
against values obtained for the polarized standard stars Velal and 
Hiltner 652; our values are within <1° from the published values, 
and the polarization percentage within 0.1 %. Table [2] provides the 
results and data source references. 



3 THE EVOLUTIONARY CODE Pegase.3 

The new version of Pegase.3 (Fioc, Rocca-Volmerange and 
Dwek, in preparation) predicts the evolution of the stellar 
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Table 2. Photometric data of the radio galaxy TN J2007— 1316 (in red on 
Fig. [5), ff(J2000) = 20 h 07 m 53 s .23, (5(72000) = -13° 16' 43.6". Comple- 
mentary data (bottom lines) are not used for fits, only given for complete- 
ness. The continuum polarization percentage accurate to within 0.1 %. 
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continuum, of the the 
quent attenuation and 
continuum by dust in 
new v ersion Pegase.3 , 
pegase. dFioc et al. 19971 ; 



metal and dust content, the conse- 
the re-emission of the UV/optical 
a self-consistent manner. With this 
u pdated from Pegase.2 (www.iap.f r/ 
iFioc and Rocca-Volmeran ge 1999b)), 
synthetic spectral energy distributions are continuously created 
from the far-UV to the sub-millime ter wavelengths at all ages 
(www.iap.fr/pegase, dFioc et al. 1997t) . We built evolutionary sce- 
narios which are consistent with the observed SEDs (and colors) 
of the full range of Hubble types at z=0. For convenience, we refer 
to these SEDs by their best fit of morphological Hubble types at 
z=0, with the understanding that they, by definition, represent the 
SEDs of galaxy representative their respective Hubble types in 
the local Universe. In Pegase.3, the chemical enrichment history 
is computed in an one-zone model, following the evolution of 
helium and most importantly, metals s uch as O, Ne, Mg , Si, S, Ca, 
C, Fe and N. We used the yields from (Marigo et al. 20 0l|) for low - 
and in termediate-mass st ars; those of (Wooslev and Weaver 19951) 
and jPortinari et al. 19981) for high-mass stars; and model W7 
of dThielemann et al. 1 986) for type la supemovae, see details 
in Fioc et al. in preparation . We assum ed that dust formed in 
circumstellar environments dDwek 1998b : in the stellar winds 
of AGB and Wolf-Rayet stars, only carbon dust is produced if 
the C/O ratio is larger than 1 and only silicate dust otherwise; in 
supernovae, as CO is unlikely to form, both carbonaceous and 
silicate grains may form but with a lower efficiency than in other 
sources of dust. We compute the dust emission by distinguishing 
between two components: star-forming HII regions and the diffuse 
interstellar medium. Monte Carlo simulations of radiative transfer 
are computed including scattering by dust grains, for geometries 
typical of disks (exponential profiles for both stars and dust) and 
elliptical galaxies (King profiles). The path lengths of photons 
before escape, absorption, or sca ttering are obtained fro m the 
method of virtual interactions dVarosi and Dwek 1999). The 
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Figure 2. An example of SED templates for the 10 1 1 M Q SO-type galaxy at 
various stages of evolution as computed with Pegase.3. The 1 fim peak of 
giant stars appears at ages > lGyr. Input parameters are: star formation rate 
= 2.10~ 3 x current gas mass perMyr, infall time scale of lOOMyrs, galactic 
winds eject gas and dust at 5 Gyr, a Kroupa IMF, and spheroidal extinction 
model (see text for details). 



attenuation is calculated for a grid of optical depths, albedos 
and asymmetry parameters, and the value at each wavelength is 
then interpolated on this grid. For each grain type and size, the 
probability distribution of temperatures is derived fo llowing the 
procedure detailed in dGuhathakurta and Draine 1989b i.e. with 
stochastic heating but assuming continuous cooling of grains. 

To model a medium of higher gas density than the classical 
interstellar medium ISM we consider the optical depth depends si- 
mul taneously on the column density NHI and the current metallic- 
ity dGuiderd oni and R occa-Volmerange 19871) . Following the for- 
malism of dDwek 1998 ) based on the respective fractions due to 
scattering K sca and to extinction A- a bs, the effective optical depth is: 



Teft = fobs + K scll ) x K.NHI 1SM .(Z C + Z si ).m H AA 

In case of denser absorbing/scattering media, the column 
density is multiplied by a AT factor where NHI KM = 6.8 x 
10 21 atoms cm -2 measured by mass unit of carbon and silicon in 
the interstellar medium of our Galaxy. 
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3.1 Pegase.3 libraries by types 

Fig.[2]shows the example of the synthetic templates for a SO-type 
galaxy at various ages as computed by the code Pegase.3. Var- 
ious synthetic libraries of SEDs evolving from ages of to 
20 Gyr are constructed to simulate the evolution of instanta- 
neous starbursts and of galaxy types of Hubble sequence. Each 
star formation scenario is defined predominantly by a set of 
four parameters: star formation law, infall timescale,fhe epoch 
of galactic winds for early types and stellar initial mass func- 
tion (see Pegase.2 read me on www.iap.fr/pegase). The param- 
eter set of each type l|Le Borgne a nd Rocca-Volmerange 20021 ; 
iFioc and Rocca-Volmerange 1999all) is chosen to predict z = 
SEDs and co lors comparable by types to local observations 
dFioc et al. 199% . The robustness of the adopted scenarios is 
checked on distant galaxies from the UV to the near-infrared 
A"-band dRocca-Volmerange et al. 2004b . They are also tested on 
faint galaxy surveys analyzed by spectral types in the mid-IR 
dRocca-Volmerange et al. 20071) and the UV -optical-nearlR do- 
mains iFioc and Rocca-Volmerange~19 99a)). Moreover they are 
used for photometric redshift prediction out high redshift (z=4) 
which shows substantial agreement with redshifts obta ined spec- 
troscopically dLe Borgne and Rocca-Volmerange 2002h . Instanta- 
neous starbursts are defined such that they form their stars over 
a period of 1 Myr. Then they passively evolve to the phase where 
their far-IR emission is dominated by the circumstellar emission 
of AGB stars and SN (the post-burst phase). The code computes 
by transfer this IR emission based on the efficiency at which the 
UV and optical photons are absorbed. Besides the SEDs by types, 
mass and age of the stellar populations are further outputs from the 
fitting procedure. Starburst libraries are computed for a variety of 
stellar initia l mass functions (IM Fs). For all scenarios, the IMF is 
taken from (Kro upa et al. 19931) after checking that other classical 
IMFs do not significantly change our predicted colours and tem- 
plate SEDs. Star formation is initiated at the redshift of formation 
that we denote as Z/ or taken as 10 when cosmic time slowly varies 
for redshift increases. While this would nominally imply an age of 
13 Gyrs (adjusted to lower redshift for Irregulars, 9 Gyrs) for local 
z=0 galaxies, this should not be taken as the average age of the stel- 
lar population as this is only when the star formation is initiated and 
not when the bulk of the stellar population forms (this depends on 
many factors of the models). Changing Zf or from z=10 to z=20 or 
30 will only vary typical ages of local templates by a small amount 
(< 0.4 Gyrs) and thus resulting in only inconsequential variations 
in the SEDs. 



3.2 Methodology 

An automatic procedure which uses a x 1 minimization searches 
for the best fit of observations. One or the sum of two spectral 
templates are from the set of synthetic libraries built for starburst 
plus a variety of SEDs fitting local Hubble type galaxies. For high 
redshift galaxies we work with template SEDs at redshift z in 
the observer's frame, by using the z=0 template SEDs reproduc- 
ing luminosities and colors of local galaxies. They are corrected 
for cosmological expansion (redshift k-correction) and for evolu- 
tion (age e-correction) using the cosmic time-z relation. For pho- 
tometry through various filters, apparent magnitudes m' A of galaxy 
type i through filter A at age t and redsh ift z are computed by 
dRocca-Volmerange and Guiderdoni 19881) : 

m\(t, z) = M\(t , 0) + (m - M)boi + k\(z, t) + e\(z, t) 



M',(/o, 0) is the local absolute magnitude by type i, (m - M\ \ 
is the distance modulus and k' A (z, t) and e' A (z, t) are the expansion 
and evolution corrections, computed from Pegase.3 synthetic tem- 
plates. 



4 RESULTS OF THE SED FITTING 

Because the two z=3.8 radio-galaxies are selected for their small 
AGN contribution to the UV/optical portion of their spectral en- 
ergy distributions as evidenced by their low level of polariza- 
tion, simple models of the thermal emission from the AGN are 
sufficient for this analysis. The mid-IR SEDs provide an up- 
per limit to the AGN component. We modeled the AGN emis- 
sion of 4C 41.17 with a black-body at T = 600 K and of 
TN J2007-1316 with the torus c mo del of dPier and Krolic 19921) 
and dKrolik and Begelman 1 988).see (Drouart et al. 20121) for de- 
tails. A more robust statistical analysis of the contribution of the 
AGN in radio galaxies will be conducted for the whole of the 
HeRGE sample. For the two radio galaxies studied here, the range 
of plausible AGN thermal continuum models has only a minor im- 
pact on the results of our stellar SED modelling and will not be 
discussed further. 



4.1 Single component model fits to the SEDs 

Fig- [Heft shows the best fits of the observed 4C 41.17 SED with 
one single burst of star formation with a duration of 1 Myr which is 
passively evolving for 25 Myr and also a single component spiral- 
like evolutionary scenario with an age of 1.6 Gyr. The starburst 
evolutionary scenario well fits the optical and far-IR (Herschel) 
data but fails to fit the region of the rest-frame 0.4-2yum (K s and 
IRAC/Spitzer) SED. The optical peak is due to stellar photospheric 
emission and the far-IR/sub-mm peak is due to a fraction of the 
Niy C photons absorbed by cold grains within environments of re- 
cent star formation. While the starbur st episode is consi stent with 
a substantial portion of the SED, dSevmour et al. 20121) for PKS 
1138-262), having to let it evolve for ~25 Myr suggests that the 
far-IR peak results from dust produced in and heated by super- 
novae and AGB stars while the optical (UV-redshifted) emission 
is due to the photospheres of recently formed mass stars. The AGN 
component, modeled as a blackbody with a temperature of 600 K, 
is compatible with data. The starburst mass, consistent with the 
best fit of the SED, is 2.7 x 10" M Q . In case of instantaneous star- 
burst, the average stellar metallicity < Z, > corresponds to the gas 
metallicity when stars were formed (null for a primitive cloud). For 
low values of < Z, > the far-IR emission is dominated by the self 
gas-enrichment of the starburst, measured by the parameter Z c +si 
which measures abundances of carbon and Silicon. The important 
point is that from the K, and IRAC/Spitzer photometry, no mat- 
ter what the exact starburst parameters are adopted, the mid-IR 
emission is not well fit. On the other hand, adopting a model with 
an older stellar population, such as an early-type spiral Sa sce- 
nario with an age of 1.6 Gyrs and having a stellar mass M, of 
3.7x10" M Q while fitting the optical, NIR and Spitzer/IRAC data 
rather well, does not fit the far-IR emission (Fig^ight). We note, 
because it is important evidence for the existence of an older pop- 
ulation in both radio galaxies, that the SED of an older population 
shows a peak at around ~ lpm which is clearly observed. However, 
despite obvious advantages to both models in fitting the SED, it is 
clear that neither single component model provides even an ade- 
quate fit to the overall SED. 
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Figure 3. The multi-wavelength SED of 4C 41.17 (dots with red error bars) compared to the sum (black line) of two components: the AGN black body at 
T t g = 600 K (dashed green line) plus a stellar component: (left) a starburst at age 25 Myr; (right) a spiral Sa at age (time since Zf or ) of 1.6 Gyr. The age, stellar 
metallicity, total stellar mass are derived from the best fitting model. Both of these single component models are missing a significant fraction of the total 
energy, either in the near-IR for the starburst or in the far-IR for the Sa-type galaxy evolutionary scenario. Vertical dotted lines mark the position of prominent 
emission lines. A second component is clearly needed to fit the SED. No intergalactic extinction ( Ay = 0) is taken into account in models. 
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Figure 4. The SED of 4C 41.17 (dots with red error bars) including the new Herschel data is best fitted (reduced x 1 = 0.96) by the sum (black line) of the 
SO-type population evolution scenario at age (time since Zf or ) of 0.7 Gyr (orange line) and a starburst with an age of 30 Myrs after the initial 1 Myr duration 
burst (blue line) from a dense (lOx standard NHIism) medium. A T e (f = 600 K black body is used to model the thermal emission from the AGN (green dashed 
line). Observations not included in stellar fits are displayed in green. As can be seen, by excluding these data, our results are relatively insensitive to the 
exact model adopted for the thermal emission from the AGN. Vertical dotted lines mark the position of prominent emission lines. No intergalactic extinction ( 
Ay = 0) is taken into account in models. 
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Figure 5. The multi-wavelength SED of TN J2007-1316 (dots with red error bars) including the most recent Herschel data is fitted (reduced ^ 2 =1. 53) by 
the sum (black line) of an old (1.2 Gyr) elliptical galaxy (orange line) plus a starburst of 35 Myrs (blue line) from a dense medium of lOx standard NHIjsm 
medium. The AGN model is the torus c model from Pier & Krolik (1992). Data points shown in green are not included in stellar fits. Vertical dotted lines mark 
the locations of prominent emission lines in the SED. No intergalactic extinction (Ay = 0) is taken into account in models. 



4.2 Two component model fits to the SEDs 

Because the one component fits are inadequate but each individu- 
ally have features which uniquely fit important regions of the over- 
all SEDs, it seems logical to attempt two component fits which re- 
alize the advantages of each single component model. To this end, 
we fit two components, one representing a recent burst of star for- 
mation and the other, a passively evolving episode of past star for- 
mation (of course constrained by the age of the universe at z~3.8). 
The best fit scenario for the radio galaxy 4C 41.17 (Fig. [4) is a sum 
of the two stellar components: a vigorous starburst observed at age 
of 30Myr from a cloud of stellar metallicity < Z, >= 5 x 1CT 4 
and mass M, = 2.2 x 10 11 M G and a second population of a mas- 
sive SO-type scenario with a short e-folding time (see Fig. [2} typ- 
ical of early-types, observed at 0.7 Gyr after its initial formation 
epoch, with a stellar mass M, = 2.0 x 10" M Q and a metallicity 
< Z, >= 5.7x 10~ 3 . Results on stellar components are weakly sensi- 
tive to the adopted form of the AGN emission (T e ff = 600K black- 
body model) in the portions of the SED which are likely dominated 
by the thermal IR emission from the AGN. 

The results on the stellar components derived from the best 
fit for TN J2007-1316 (Fig. [5j are generally similar to those for 
4C 41.17, namely there is the need for a young starburst and an 
older stellar population similar to that of a local massive galaxy 
evolved by a redshift of 3.8. Specifically, the fit requires a starburst 
of very low metallicity evolved for 35 Myrs with a total stellar mass 
of M, = 8. 1 x 10 10 M Q and a second component consistent with that 
of an evolved massive early- type galaxy observed at age of 1.2 Gyr 



(after z for ) with the average stellar metallicity (< Z, >= 1.5 x 10 2 ) 
and a significant stellar mass (M* = 9.4 x 10" M ). 

4.2.1 A massive post-burst in distant radio galaxies 

For the two z = 3.8 radio galaxies, fitting simultaneously the op- 
tical 1500A rest-frame) and crucially the far-IR SED provided 
by the Herschel/submm data suggests that both galaxies are under- 
went a massive starburst several tens of Myr ago ( a post-burst). The 
modeling of the far-IR emission depends on the metal-enrichment 
of Carbon and Silicon, Zc+si, while the observed UV-optical emis- 
sion is dependent on the amount of absorbed photospheric emis- 
sion intrinsic to the young stellar population. The main character- 
istics of the fitted young post-burst scenario are summarized in Ta- 
ble. [3] The best fit is a result of having to balance the far-IR and 
optical luminosities and spectral shapes, measured by the reduced 
X' rei minimum. A younger age increases the number of Lyman con- 
tinuum photons but decreases the metal-enrichment, favoring the 
direct photospheric emission relative to dust emission. Older ages 
have the opposite effect. Because of this dependence, the starburst 
mass and age are well constrained and unique. All other parameters 
of Table. [3] (Z c+Si metallicity , supernovae number, high density 
value) are not free parameters, but derived from chemical evolu- 
tion. We recall that the mean metallicity of stars < Z* >, as indi- 
cated in the legends, only traces the gas metallicity when stars were 
formed. In case of starbursts, this parameter is the initial gas metal- 
licity and remains constant with age while gas metallicity evolves 
and is largely dominated at age of 30 Myrs by the intrinsic gas- 
enrichment due to the passive evolution of the starburst scenario. 
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Table 3. Characteristics of the young post-burst at z=3.8: age, stellar mass M*, initial star formation rate, bolometric luminosity htoi, dust/bolometric lumi- 
nosity ratio, carbon+silicon Zc+Si mass fraction, column density factor K (NHI)= K X NHIvsm), and number of type-II SNe. 
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Table 4. Characteristics of the old-star population at z = 3.8: age, stellar mass M„, current star formation rate, bolometric luminosity, dust/bolometric 
luminosity ratio, Carbon+Silicate dust mass fractions, SNIa number, SNII number. 
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In the radio galaxy 4C 41.17 (Fig. [4] and Table. [3} the new 
Herschel/submm observations are well modelled by the cold grain 
emission of the highly massive and metal-enriched (Z c+- s, =0.026) 
stellar environments. This post-burst population was triggered 
30 Myrs earlier in a starburst with a duration of 1 Myr (instanta- 
neous) with a star formation rate of 2.6 x 10 s Af G .yr~' (current stel- 
lar mass corrected for already died stars at age of 30 Myrs). The 
rapid evolution of most massive stars (a few 10 7 yrs) means that 
they have already evolved to SN and AGB stars when the post- 
burst is observed, enriching the interstellar medium. This implies 
that supernovae explosions (usnii = 7.1 x 10 7 by Myr) and AGB 
stars with circumstellar envelops are dominant. A dense medium 
with the high column density (a factor K = lOx the column density 
measured in the Galaxy interstellar medium NHI /-SM ) is required by 
the far-IR data, implying 10 22 ~ 23 atom s cm -2 . Such column densi- 
ties have been found in radio galaxies (Mu llanev et al. 2010). This 
post-burst in the far-IR dominates the total luminosity (96 %). In the 
radio galaxy TN J2007- 1316 (Fig.[5]and Table. [3}, the starburst age 
is 35 Myrs with a mass of 8.1 x 10'° M Q and an initial metallicity of 
zero, likely a primitive cloud. This suggests that most of the neces- 
sary metal enrichment to explain the far-IR emission (Z c+S ,=0.025) 
was due to self-enrichment of the young stellar population. 



4.2.2 The importance of the 1 pm peak 

Fig. [4] and Fig. [5] show the typical energy distribution of an 
old stellar population at wavelengths of the A"-band and through 
SpitzerflRAC filters (from the 4000 A discontinuity to the typi- 
cal 1 pm peak in rest-frame). The evolution of early-type galax- 
ies are characterized by a short (^1 Gyr) intense star formation 
episodes. This old population is seen in both targets perhaps sug- 
gesting it is a generic feature of radio galaxies. For both radio 
galaxies, this component has an age that is a significant fraction 
of a Hubble time at z=3.8. We recall that we have characterized 
the ages as the time from a formation redshift Zf or , which for 
4C 41.17 and TN J2007-1316 is roughly 0.7 Gyr and 1.2 Gyr re- 
spectively. However, acceptable fits allow for a range of age at least 
a factor of 2 around these nominal values. We also note that the 
bulk of the stars in this old population formed more recently than 
the values listed. This peak in the SED arising from evolved gi- 
ant stars only evolves slowly and has been identified at all red- 
shifts. This result supports our analysis of the Hubble A"-band 



J resttam, (A) 



10"" 




. ^ ' 1 ^ g 1 1 1 1 1 1 11 


10"' 8 






;> io" 19 

i 

\ 






5 10" zo 






iq -8 ' 






io~ S2 


TN J2Q07--316 
vL = 1.7 J , rcdshitt = J.B40 ; Av=0.0D 

■ E2 ; QC(?1= 1400 tvivrs ; l/- = 2.bo+ '■ 1 !,■„ : <ZO = 0.0257 
I ; oc.p.?= 40 h,\rs ; M* = 3.1e+10 M 3 : < Z* > = 0.0 ' 00 

■ AGN torus model, i = 3!deg 



10 4 10 5 10 6 10 7 

» (a; 



Figure 6. Best fit (reduced ^ 2 =0.98) of the multi-wavelength SED of 
TN J2007-1316 (dots with red error bars) with the sum (black line) of a 
40 Myrs starburst blue line) plus an early -type galaxy (1.4 Gyr) after an 
episode of galactic winds (orange line) ejecting gas and dust, they stop any 
further star formation. Masses and ages vary by ^ 50% compared to results 
of Fig. m 

diagram dRocca-Volmerange et al. 20041) . predicting already mas- 
sive and old galaxies at z > 4. Table. [4] presents the characteristics 
of the evolved old component. Stellar masses are respectively 2x 
10" M for 4C 41.17 and a 10 12 M o for TN J2007-1316. These 
masses are higher than that of the young starburst mass and per- 
haps suggest previous episodes of merging at earlier epochs. The 
contribution of the AGN to MlPS/Spitzer data, modeled by a sim- 
ple black-body law at T e ff = 600 A* as shown on Fig. [4] (dashed 
green line), and by a more sop histicated torus c-m odel (adapted 
from dPier and Krolic 19921) by dDrouart et al. 20121) ) for the radio 
galaxy TN J2007-1316 (Fig.[5]dashed green line), does not change 
our conclusions: for these two radio galaxies, the AGN emission is 
not significant in the near-infrared domain at z = 3.8. 



5 DISCUSSION 

The SEDs of the radiogalaxies 4C 41.17 and TN J2007-1316 
are revealed with two stellar, respectively young and old compo- 
nents and a simple AGN model. The spectral synthesis used to fit 



the data is based on a variety of template libraries including in- 
stantaneous starbursts and evolutionary scenarios designed to fit 
the SEDs of local Hubble-type galaxies. Each scenario is charac- 
terized by a star formation time-scale: ~1 Myr for the instanta- 
neous starburst, ~ 1 Gyr for elli pticals to ~10 Gyrs for spirals (see 
dRocca-Volmerange et al. 20041) . In case of gas-dependent SF laws, 
star formation time-scales are regulated by the accretion and galaxy 
scale outflows of gas. One template library is built for each sce- 
nario. For starbursts, various standard IMFs are modeled in these 
template libraries. Because in evolutionary models as Pegase.3, 
IMFs are all normalized to 1 M Q and is held fixed throughout the 
evolution, changing the IMF does have an effect on relative rate of 
metal-enrichment and dust fraction, but has only a minor effect on 
the best fit results of the SEP. Although onl y weakly constrained, 
all of our best fits favor a dKroupa et al. 19931) initial mass function. 
The most robust determinations from our fitting analysis are the 
relative ages and total stellar mass of each best fit scenario. For the 
two galaxies, our results suggest that shape of the SEDs is due to a 
vigorous recent starburst superimposed on a relatively evolved pop- 
ulation whose formation was initiated at early epochs (zf or > 10; 
but this is not when the bulk of the stellar population was formed, 
only when the star formation was initiated within the context of the 
scenarios used to model the SEDs). For the young starburst com- 
ponent, our results suggest an age of ~30 Myr with a stellar mass 
of =i 10" M . The two peaks respectively seen in the optical and in 
the fazIR/Herschel are the data which constrain this fitted compo- 
nent the most. This is because of the timescale for dust production 
in the postburst AGB and SN populations. Of much lower precision 
is the age of the older component which is mainly constrained by 
the peak at rest-frame 1 /jm due to slowly evolving giant stars. This 
slow evolution means that a wide variety of ages could in princi- 
ple fit the data roughly equally well, but all best fit models suggest 
relatively old ages (100s of Myr). The best fits are that of an SO 
galaxy with an age of 700 Myr for 4C 41. 17 and an elliptical galaxy 
scenario with an age of 1200 Myrs for TNI 2007-1316. Both are 
early-types with a ~1 Gyr time-scale and masses of = 10' 1-12 M B . 
To check the sensitivity of the 4C 41.17 results to ages of the mod- 
els, we tested the relative variations in x 2 for different models. We 
maintain an instantaneous starburst as first component. Replacing 
the second component by another starburst of 1 Myr duration in- 
duces a x 1 variation of 10% and an age reduction of 80% with 
deeper balmer lines. Replacing the second component by an later- 
type (Sa) scenario with star formation scale of a few Gyrs induces 
ax 2 variation of 5% and an age increase of 120%. 

For TN J2007-1316, no other scenario, excepted ellipticals, 
is acceptable within a^ 2 variation of 400%. We also test (Fig. [6} a 
second component of an early-type model suffering galactic winds 
at early epochs (at age of 1 Gyr). The 1/rni peak is then better fitted 
to IRAC data at age of 1 .4 Gyr and summed with a younger (40 
Gyrs) starburst. For the two populations, ages do not vary more 
than 50%. 

In summary for these two powerful radio galaxies, the young 
post-burst and an evolved component of a short time-scale (~1 Gyr) 
are the best solutions. The timescale for both the increase in the 
metallicity with time and relative fraction of dust in the ISM which 
produced the far-IR emission are calculated in a self-consistent 
manner based on stellar yields and the production of dust by stars. 
The relatively poor time resolution in fitting the broad band contin- 
uum with stellar synthesis models for high redshift galaxies does 
not allow us to discuss anything but simple evolutionary models 
(not for example successive short bursts which would look more 
like continuous star formation in the models). Long star-formation 
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Table 5. Starburst/old population ratios 



Radiogalaxy ages stellar mass SNII number factor K 

4C41.17 30:700 2.2:2.0 71:1.9 10:1 

TNJ 2007-1316 35:1200 8:94 30:6.7 10:1 



time-scales results in two cumulated effects: increasing metal and 
dust fraction and simultaneously decreasing the photospheric emis- 
sion of massive stars. Because of these effects, the parameter-range 
which simultaneously fits both the rest-frame optical and infrared 
portion of the SEDs is relatively narrow. That is why it is important 
to do this analysis on multiple sources, which if give consistent re- 
sults, means that the general characteristics of the best fits are likely 
to be robust. 

The general conclusions we can make are that the masses are 
likely in the range 10' ~ u M o for recent starbursts and 10"~ 12 M o 
for the early-type evolutionary scenarios used to fit the data. Star- 
bursts and evolved galaxies, have mass ratios of 1:1 to 1:10 typical 
of major mergers, so that any mass interchange will not signifi- 
cantly affect results. Moreover all mergers are gas-rich inducing an 
efficient star formation. 

5.1 Star formation process in 4C41.17 

What is triggering so intense star formation in distant radio galax- 
ies? One of the most striking feature of these distant SEDs is the 
double peak of the starburst activity, respectively in the far-IR Her- 
schel from dust and the optical, emitted from star photospheres. 
The instantaneous starburst model evolves to the age of ~30 Myrs 
old likely dominated by supernovae and AGB stars. The huge stel- 
lar mass =! 2.210"M o corresponds to a rapid metal-enrichment 
from a no-primitive cloud (< Z, >= 5 10~ 4 ). The emission is is- 
sued from a gas of high density, characterized by the factor K = 
10, describing a ten time denser ISM that the classical interstel- 
lar medium. A model has been already pro posed for 4C 41.17, a 
template of the s o-called alignment effect (Chambers et al. 19901 ; 
iMilev et al. 1992]) . Based on the jet-cloud interaction inducing star 
formation by sh ocks, suggested by m any authors from 1980, it 
was revisited by jBicknell et al. 2000|) in the light of HST and ra - 
dio emission line data Jvan Breugel et al. 19981 ; iDev et al. 19971) . 
The authors conclude to a set of parameters: baryonic mass of a 
8xl0 10 Mq, a few Myrs star formation time-scale, a dynamic radio 
galaxy age of = 30 Myrs and the high density of ^ 1-10 cirr 3 . All 
these values are compatible with our stellar population model, fa- 
voring the jet cloud interaction model for triggering star formation 
at high redshift. Another favorable argument for this scenario, at a 
much lower rate, is the star formation by jet-cloud interaction dis- 
covered in the core of NGC 1068 w ith the high angular resolution 
of NaCo/VLr texposito et al. 201 ll) . 

The dense medium measured by a column density factor 
K=\0 meaning 10 22 ~ 23 atoms cirr 2 only seen in the starburst com- 
ponent requires justifications. A better spatial resolution associated 
to a more refined 3D spectroscopy with the new generation of in- 
struments (MUSE/VLT) will solve the location of star formation 
activity and the feedback measurement. 

The coincidence of the short timescales of both starburst and 
jet-interaction could only be due to a similar causal origin. In that 
case, the post-burst we are seeing at works is only an episode 
of an active merger process. Many other arguments would fa- 
vor massive gas-rich mergers at high redshifts. The discovery 
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of huge CO comp onents toe Breuck et al. 20051 ; lingel et al. 2o"l(il : 
llvison et al. 20121) already suggest evidences of major mergers in 
such powerful distant radio galaxies. Another convincing remark 
in favor of gas-rich mergers is that to form so huge stellar masses 
in short episodes requires large volumes of high gas-density. Ta- 
ble. |5]gives the characteristics of the intervening masses, with mass 
ratios typical of major mergers of young and evolved components 
of various metallicity, occasionally of primitive gas, implying sev- 
eral star generations through merger processes. 

5.2 The old early-type stellar population 

Local e arly-type galaxies a re generally not luminous in the 
far-IR dXilouris et al. 20041) . which has been interpreted as 
the result of galaxy win ds driven by the evolved populations 
dMafhews and Baker 19871) The other possible mechanism for driv- 
ing winds is the mechanical and radiative energy from intense 
starbursts. For example, the local archetypal starburst galaxy 
M82, galactic winds were recently studied with HerschelfPACS 
dContursi et al. 20121) . The authors propose that cold clouds en- 
trained in the galaxy disk outflow. At high redshifts, star formation 
activity is intense enough at the earliest epochs to induce intense 
winds which may also entrain the molecular and neutral atomic hy- 
drogen driving it out of the galaxy. This assumption will possibly 
be statistically analyzed on the totality of the HeRGE sample. 

The new important result is the emergence in the mid-IR (op- 
tical rest-frame) of an (already evolved) stellar component. Seen 
through the specific window (IRAC filters). The age and huge stel- 
lar mass of this old population is calibrated on the SEDs of two 
selected radio galaxies, in particular when galactic winds are im- 
posed. At age of ~ 1 Gyr, the evolved stellar population is giant 
stars building the well-known peak at lyum (rest-frame; Fig.|2j- All 
these results favor the presence of highly massive early type galax- 
ies initiated at early epochs Zf or >10, which must have grown very 
rapidly. Because the insufficient time resolution, the declining con- 
tinuous star formation laws, as used in Pegase.3 and other codes 
could also be considered as a series of successive starbursts with 
intensities declining with time. From this point of view, which is 
also consistent with the data, the hierarchical assembly history of 
most ma ssive (~8xl0 12 M Q ) haloes in a ~3 Gpc 3 volume as pro- 
posed in dLi et al. 2007b to form populations of z=6 quasars would 
be compatible with our results. 



5.3 The link with radio-quiet massive galaxies 

The discovery of these z = 3.8 old stellar component in power- 
ful radio galaxies with a superposed starburst in a high density 
medium, warrants a comparison with radio-quiet early-type high- 
z galaxies. These radio quiet galaxies have med ian stellar masses 
of ~4x 10 11 M at z=2.0-2.7 lkriek et al. 20o3) . Their mass esti- 
mates are likely underestimated because observations were through 
the Gemini K, filter, missing the population of dominant rest-frame 
K-band emission redshifted in the Spitzer filters. However the com- 
parison with powerful radio galaxies is justified because their re- 
spective masses are of the same order. Moreover deep and high- 
resolut ion images obtained wi t h RST/NICMOS and /NIC2 and 
Keck dvan D okkum e t al. 20081 ; IStockton et al. 20 08) are associ- 
ated to the presence of massive disks of old stars at high redshift of 
remarkably small sizes with a median effective radius, r e ff=0.9kpc. 
Cosmological surface brightness dimming has a strong impact on 
the perceived morphology of distant galaxies and likely explains 



generally why disks are basically undetectable at z> 1.5. The decre- 
ment in surface brightness due to cosmological dimming is of 7 
magnitudes at redshift 3.8. Only the central disk and spheroidal 
components are of sufficient ly high surface brig htness to be de- 
tectable even in deep images dStockton et al. 20081) . 

If radio-load galaxies are a subpopulation of the massive 
radio-quiet galaxies which are undergoing an episode of vigorous 
star formation, it may well be that the powerful AGN illuminate the 
surrounding gas and dust in the disk over larger scales than would 
not otherwise be visible due to surface brightness dimming. Even 
the continuum emission would be affected due to scattering of the 
AGN continuum by dust. Moreover, in the population of radio-quiet 
massive galaxies, the stellar population does not have the massive 
young stellar component we have observed in IR-luminous pow- 
erful radio galaxies, and therefore would be less luminous in both 
the far-IR and in the optical perhaps below the threshold of the 
detectability. The massive evolved starburst in IR-luminous pow- 
erful radio galaxies contribute to the illumination of the disk over 
larger scales, sticking out of the surface brightness threshold. Ob- 
serving extended emission preferentially in radio galaxies confirms 
the coexistence of the radio and starburst activity in the short star 
forming episode perhaps initiated by gas-rich mergers. Of course, 
the jet would also have an impact through its interaction with the 
surrounding gas, compressing and heating it. This shaping of the 
gas by the radio jet may also lead to stronger extended emission in 
both the emission lines and continuum. 

More knowledge on kinematics and apparent morphology 
strongly modified by distance effects are also required for justify- 
ing the spheroidal structure of early-type galaxies. In case of radio 
galaxies, all these properties have to be associated to the presence 
of a supermassive black hole. Future studies will focus on thelumi- 
nosity relationship between AGN and star formation activity of a 
larger sample of HeRGE galaxies. 



6 CONCLUSIONS 

From the Pegase.3 evolutionary spectral synthesis of two z = 3.8 
radio galaxies, the SED fitting on a large wavelength coverage 
(UV/optical to far-IR/submm), identifies two distinct stellar com- 
ponents. The originality and specificity of the analysis is to work in 
the observer's frame by using robust local templates, corrected at 
high z for the Universe expansion (k-corrected) and evolution due 
to distance (e-corrected), following scenario librairies by types of 
the code Pegase.3. At any time, synthetic templates are continu- 
ously computed with coherent UV/optical stellar emission and dust 
absorption reemitted in the far-IR. Instantaneous (lMyr) starbursts 
and a variety of scenarios of the Hubble sequence, varying with 
1-10 Gyrs star formation time-scales, are considered for building 
template libraries. Best SED fits are derived in the observer's frame 
by &x 2 algorithmic procedure on the sum of two stellar components 
plus a simple AGN model. 

Main results are similar for the two galaxies at z = 3.8: 

• one single component is unable to explain the complete 
spectral energy distribution, whatever types and ages. 

• the sum of two stellar components give the best SED fits, 
following: 

• A 1 Myr starburst, observed at age of 30 Myr, forming a huge 
stellar mass = 10" Mo in a dense (lOx standard ISM) medium, 
dominant in the far-IR Herschel to submillimetric and in the optical 
domains. 

• The discovery in the K-band to Spitzer/mid-IR domains 
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(lyum rest frame) of an already evolved (=d Gyr at z = 4) popu- 
lation. 

• The evolved stellar component is compatible with an early- 
type galaxy ( SO-Elliptical) initiated at early epochs (z/ or > 10), 
with possible galactic winds. 

• The mass of the evolved component is huge of 10"~ 12 M G 
for the two galaxies at z —4. 

• The possibility of a gas-rich merger (joined to a possible 
jet-cloud interaction) is favored. 

Over the longer term, we intend to include a more refined 
modeling of the AGN contribution to the overall SED in the Pe- 
gase.3 evolutionary synthesis models and to model the optical, Her- 
schel/submm and Spitzer SED data of the HeRGE sample over its 
entire redshift range (1 < z < 5). 
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